function [J] = Q4_J(r,s,xy)
Nrs=0.25*[-(1-s) 1-s 1+s -(1+s);
         -(1-r) -(1+r) 1+r 1-r];
J = Nrs * xy;
end